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O ABSTRACT 
H 

The relative roles of metallicity and surface gravity on the near-infrared spectra of late-T 
d brown dwarfs are not yet fully understood, and evolutionary models still need to be calibrated in 

order to provide accurate estimates of brown dwarf physical parameters from measured spectra. 
I The T-type brown dwarfs Eps Indi Ba and Bb forming the tightly bound binary Eps Indi B, which 

^ orbits the K4V star Eps Indi A, are nowadays the only such benchmark T dwarfs for which all 

important physical parameters such as metallicity, age and mass are (or soon will be) known, 
as We present spatially resolved VLT/NACO images and low resolution spectra of Eps Indi B in 

the J, H and K near-infrared bands. The spectral types of Eps Indi Ba and Bb are determined 
by direct comparison of the flux-calibrated JHK spectra with T dwarf standard template spectra 
and also by NIR spectral indices. Eps Indi Bb is confirmed as a T6 while the spectral type of 
Eps Indi Ba is T1.5 so somewhat later than the previously reported Tl. Constrained values for 
surface gravity and effective temperature are derived by comparison with model spectra. The 
evolutionary models predict masses around ^53 Mj for Eps Indi Ba and ~34 Mj for Eps Indi Bb, 
. slightly higher than previously reported values. The suppressed J-band and enhanced K-band 

flux of Eps Indi Ba indicates that a noticeable cloud layer is still present in a T1.5 dwarf while 
^ no clouds are needed to model the spectrum of Eps Indi Bb. 
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1. Introduction 

Multiple stellar systems with one component 
being a main sequence star and at least one other 
component having a mass below the hydrogen- 
burning limit, and hence residing in the domain of 
brown dwarfs are the key to calibrate and bench- 
mark brown dwarf evolution models. Since all 
components are expected to be coeval and of the 



same chemical composition, age and metallicity of 

the main sequence primary determined by well es- 
tablished methods also apply to the brown dwarf 
companion(s). Furthermore, sufficiently tight bi- 
nary systems would allow one to directly measure 
the system mass by determining the orbit and even 
to measure the individual component masses via 
monitoring of radial velocities. Hence, the mea- 
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sured spectral energy distribution (SED) of such 
a companion can be used to calibrate the effective 
temperature and surface gravity of evolutionary 
models and provide the required benchmarks. 
One such precious object, Eps Indi B, was re- 



solved into a binary T1/T6 dwarf (McCaughrean 
et al.|[2004 l orbiting the K4.5V primary Eps Indi 
A during the commissioning of the Simultaneous 



Differential Imager (SDI, [Close et al., 2005} at 
the VLT. Initial mass estimates based on spec- 
tral properties and an assumed age of 1.3 Gyr were 
47±10Mjup and 28±7Mjup for Eps Ind B a and Bb, 
respectively ( McCaughrean et al. 2004 ) . Another 



T6 dwarf companion, SCR 1845 B, was recently 



found by Biller et al. (20061 around the MS. 5 pri 



mary SCR 1845 A using the same instrument and 
characterized by NIR imaging and spectroscopy 
by [Kasper et aT] (|2007|). Both, Eps Indi B and 



SCR 1845 B, are extremely nearby (closer than 4 
pc), so near that they can efficiently be observed 
with modern Adaptive Optics instruments despite 
their intrinsic faintness and rather small orbital 
distance from the primary. Other benchmark ob- 



2002) 



jects include HN Peg (Leggett et al. 20081, Gl 



570D (Saumon et al. 20061 and HD 130948 BC 



(Dupuy et al. 2008 Potter et al. 2002 Goto et al 



It is now widely accepted that, besides Tog 
(a function of bolometric luminosity and radius 
which are both a function of age), surface gravity 
(a function of mass and radius) and metallicity 
are centrally important in determining the spec- 



tral energy distribution of late-T dwarfs ( Saumon 
et al.||2007| |Liu et al.|[2007| . The effects of metal- 
licity and surface gravity on the NIR spectra are, 



however, not yet fully understood ( Burgasser et al. 
2006a Liu et aL]|2007 ). Calibration of theoretical 
models against observations has started with the 



work by Burgasser et al. (2006a I; Saumon et al. 



(2006); Leggett et al. (2007). Both mass and age. 



i.e., the main underlying physical parameters for 
any modelling of brown dwarf evolution and at- 
mospheres, are, however, still highly uncertain for 
most known T dwarfs. 

Precise mass estimates can only be derived for 
brown dwarf, which are members of binary or mul- 
tiple systems, and for which orbital parameters 
have been determined. Currently, Eps Indi Bb is 
the only late T dwarf for which physical parame- 
ters, and in particular its mass, will soon be known 



with high precision. While the Eps Indi B binary 
is separated from the Eps Indi A by 1500 AU (sev- 
eral arcminutes on sky at the distance of 3.626 pc), 
the two brown dwarfs are separated only by about 
0.7" corresponding to a projected spatial separa- 
tion of 2.6 AU (iMcCaughrean et al.||2004|). This 



suggests that the system mass might be measured 
from orbital motion within a few years. From mea- 
surements on Eps Indi A, the system age is rea- 
sonably well determined to be in the range 0.8 to 
2 Gyr ( Lachaume et al.|1999 ) , and a near-to-solar 
metallicity ([Mg/H] = - 0.05,[A1/H] ^ -0.04 [Fe/H] 
= -0.06) is reported by Beirao et al. (2005). 



2. Observations and Data Reduction 

Eps Indi B was observed with NACO mounted 
on UT4 of the VLT on 23 May 2008 (imaging) and 
9 June 2008 (spectroscopy). The Adaptive Optics 
used 90% of the light from the binary as a guide 
star for its infrared wavefront sensor while 10% 
of the light were directed towards the scientific 
camera. In this way, we obtained well corrected 
images that easily resolve the 0.58" binary (see 
Figure[l]) . 

2.1. Imaging 

For the imaging observations, the S13 camera of 
NACO was used in the J, H- and Ks near-infrared 
filters, with a pixel scale of 13.20 ± 0.05 mas per 
pixel. The target was observed using a standard 
dithering procedure integrating for three minutes 
on source in each of the bands. 

Image data reduction followed standard proce- 
dures including sky removal, bad pixel cleaning, 
and flat fielding. Differential photometry and as- 
trometry were done on the reduced images by si- 
multaneously fitting a model PSF profile to the 
binary. The individual magnitudes for Eps Indi 




Fig. 1. — NACO broad band JHK images of the 
0'.'583 binary Eps Indi Bb. North is up. East is 
left. 
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Ba and Bb were finally calculated from the rela- 
tive magnitudes and the integrated 2MASS magni- 
tudes of the unresolved system (J=11.91, H=11.31 
and Ks=11.21). The error introduced by the dif- 
ference between NACO and 2MASS instrument 
responses has been calculated from the spectra, 
as described in the next section, and has been 
taken into account. Tablejl] displays 2MASS pho- 
tometry of Eps Indi Ba and Bb. On 23 May 
2008, the separation between the two components 
was 0'.'583±0'.'003, and the position angle was 
153?9±0?5. The errors of the astrometric posi- 
tions and the flux are the 95% confidence levels 
provided by the MATLAB fitting algorithm, and 
the error on position angle is introduced by the 
uncertainty of the NACO FoV orientation. 

2.2. Spectroscopy 

NACO was used in long-slit spectroscopic mode 
to obtain R « 400 spectra in the J and HK spec- 
troscopic filters (modes S54-4-SJ and S54-4-SHK, 
2nm/px pixel dispersion, S54 camera with 54.3 
mas/pixel). By turning the Nasmyth adaptor ro- 
tator, the slit was aligned with the binary axis. 
Dithering the binary at two positions along the 
slit, exposures of 26 minutes total in each of the 
modes were obtained. A telluric calibrator (HIP 
113982, spectral type AOV) has been observed in 
a similar manner just after the object at a simi- 
lar airmass. The effects of differential slit losses 
between Eps Indi B and the spectral calibrator 



(Goto et al. 20031 were mitigated by observing 
both properly centered on the slit with about the 
same orientation of the parallactic angle with re- 
spect to the dispersion direction. A spectrum of 
the NACO internal Tungsten lamp was recorded 
for wavelength calibration. 

Spectroscopic data reduction followed standard 
procedures. After sky removal, bad pixel clean- 
ing and fiat fielding, the spectra from Eps Indi 
Ba and Bb were extracted and wavelength cal- 



Table 1: 2MASS photometry of Eps Indi Ba and 
Bb from 23 May 2008 (JD 2454609.9). 



Star 


J 


H 


Ks 


Ba 


12.29±0.03 


11.52±0.03 


11.36±0.03 


Bb 


13.23±0.03 


13.21±0.03 


13.47±0.03 



ibrated. The halo of Eps Indi Ba at the po- 
sition of Bb was estimated at its opposite po- 
sition with respect to Eps Indi Ba. The spec- 
tra were divided by the telluric calibrator and 
multiplied by a Kurucz template spectrum of 
Vega (http:/ /kurucz. harvard.edu/stars/VEGA/) 
smoothed to the actual resolution of the data in 
order to remove the hydrogen absorption features 
of the AOV telluric calibrator. Synthetic magni- 
tudes were calculated as described in |Kasper et al.| 
(2007) in order to determine and apply the cor- 



rection factors needed to convert the fiux ratios 
from the NACO to the 2MASS photometric sys- 
tem. Flux calibration of the spectra was finally 
established using the calculated 2MASS magni- 
tudes of both components (see tabl e[l]) and the 
2MASS zero magnitude in-band fiux (Cohen et al. 



2003D . 

It was verified that the errors introduced by 
the difference of the isophotal wavelengths for Eps 
Indi B and and the spectral calibrator do not ex- 
ceed the photometric errors. We also estimated 
noise from the standard deviation of the individ- 
ual spectra that were generated during the data 
analysis and averaged to obtain the final result. 
The SNR for each individual data point is around 
30 for Eps Indi Ba and most parts of Eps Indi Bb 
spectra. In the K-band and in the depression be- 
tween the J-band peaks, the SNR of the Eps Indi 
Bb spectrum drops to 3-10. The deduced NIR 
spectrum is consistent with the optical to mid-IR 



spectra presented by King et al. (20081. 



Results 



3.1. Spectral Types 



The most direct way to determine the spectral 
type of a star is direct comparison to spectral stan- 
dards. Figures[2] and [3] plot the NIR spectrum of 
Eps Indi Ba and Bb (black lines) on top of T dwarf 



et al. 



spectra obtained by Knapp et al. (20041; Strauss 



( [T999| ); |Leggett et al.| ( |2000[ ); |Chiu et al 



(2006) and defined as spectral standards by Bur 



Igasser et al. ( 2006b I . The NIR spectrum of Eps 
indi Ba fits somewhat in between the standard Tl 
and T2 spectra. Using the spectral indices defined 



by Burg asser et al.| ( |2006b| ), three (H2O-J, H2O-H, 
CH4-K) suggest Tl while CH4-H suggests T2 and 
CH4-J even T3. Hence, we slightly refine the spec- 
tral type of Eps Indi Ba to T1.5. The spectrum of 
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Eps Indi Bb fits the T6 standard well. 



3.2. Fitting Models to the Eps Indi Ba+b 
Data 

Assuming tlie Eps Indi B components iiave tlie 
same metallicity (near solar, Beirao et al.' ('2005')) 
and age (~ 0.8 — 2 Gyr, Lachaume et al. (1999)) 



allows us to explore the accuracy of basic brown 
dwarf evolution and atmosphere theory. Evolu- 
tionary theory provides the mapping between, on 
the one hand, mass, age, and radius and, on the 
other, Tgff and gravity (g) (for a given metallic- 
ity). For given Teff and gravity, we have a pre- 
dicted radius, and with the known distance we 
can then make absolute comparisons between the 
meaured spectra and theoretical models. For the 



evolutionary theory, we use Burrows et al. ( 1997 1 



For the spectral and atmosphere theory, we use the 



approach described in Burrows et al. (20061 and 



Sharp fc Burrows (20071. The goal is to deter 



mine Teff and gravity for each object, and then to 
determine whether a consistent theoretical picture 
emerges. In particular, we can ascertain whether 
the models are consistent with the assumption of 
coevality. 

These current models have some limitations. 
First, since Eps Indi Ba is clearly near the L/T 
transition region above a Teff of ~1100 K, it 
is necessary to incorporate the effects of sihcate 



clouds ( 


Ackerman & Marley 


2001 


Tsuji 


2002 


Burrows et al. 


2006 


). The modeling of clouds 



(e.g., their modal particle sizes, spatial extent, op- 
tical properties, and degree of patchiness) is quite 
problematic. Nevertheless, we use the "Case E" 



cloud model described in (Burrows et al. 20061 
with 15 /im spherical forsterite grains of which the 
cloud is assumed to be comprised. We initially 
also tested 10 and 20 /zm particles but did not 
achieve better fits. Second, though the absorp- 
tion feature due to the hot bands of methane from 
^^1.6 to ~1.75 /im is a defining signature of brown 
dwarfs, the strengths of these bands are not well- 
determined theoretically, and have not been mea- 
sured in the laboratory. Since it is as good as any 
credible formulation described in the extant lit- 
erature and though it clearly underestimates the 
strength of these bands, we use the prescription 
described in Sharp & Burrows (20071. Finally, 



the wings of the K I resonance absorption doublet 
at ^0.77 fxrn determine the continuum of brown 




Fig. 2. — Spectrum of Eps Indi Ba (black line) 
plotted on top of Tl-3 dwarf spectral standards 
obtained from Sandy Leggett's L and T dwarf data 
archive. 
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X[|im] 



Fig. 3. — Spectrum of Eps Indi Bb (black line) 
plotted on top of T5-7 dwarf spectral standards 
obtained from Sandy Leggett's L and T dwarf data 
archive. 
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dwarf spectra out to as far as ^-^1.05 fim (Burrows 
et al.|[2000 l. The shape of these wings, in particu- 
lar the effective cutoff wavelength, is still a subject 
of active investigation. 

Given these caveats, we have identified the 
range of spectral models which approximately re- 
produce the observations. For Eps Indi Ba, we 
ran grids of models spanning the Toff range 1150- 
1300 K with a resolution of 50 K and the log(g [cm 
s^^]) range 4.9-5.3 with a resolution of 0.1. For 
Eps Indi Bb, we ran grids of models spanning the 
Teff range 850-950 K with a resolution of 25 K and 
the \og{g) range 4.7-5.1 with a resolution of 0.1. 
The models were than compared to the measured 
spectra by visual inspection looking for the small- 
est residuals between data and model over the 
complete spectral region. We further calculated 
least-squares residuals between data and model 
and found that the results agree very well with 
the by-eye fits. We did not attempt to use more 
sophisticated numerical methods such as statis- 



tics ( [Gushing et al.|2008[[&]rgasser et al.|2007[ ) to 

quantify the goodness-of-fit because current brown 
dwarf spectral theories have yet too many unquan- 
tified systematic problems that make a rigorous 
numerical treatment difficult and to some extent 
still arbitrary. Masses and ages were derived for 
the three best fitting models from the evolutionary 



tracks of Burrows et al. ( 1997 1 by linear interpola- 
tion. We then determined the final mass and age 
estimate by calculating a weighted average of the 
three best fitting models assigning weights of 1, 0.5 
and 0.33 in order of preference. Errors on these fi- 
nal estimates were derived from the mass and age 
variations introduced by our log(.g) step size of 0.1 
which dominate over the variations introduced by 
the chosen Tcs resolution. 

For Eps Indi Bb, no cloud models were nec- 
essary and we derive well matching fits for (Tog 
/ log{g)) of (900 / 5.0), (925 / 5.1) and (875 / 
4.9) in order of preference. The best-fit values of 
Teff and g are correlated, with the higher gravities 
paired with the higher TcffS. Using the evolution- 
ary tracks, the corresponding masses and ages are 
(34/,Mj / 1.5 Gyr), (38Mj / 1.8Gyr) and (30Mj / 
1.1 Gyr). Using the weighted average described in 
the previous paragraph, we finally estimate mass 
and age of Eps Indi Bb of 34±5 Mj and 1.5^°-^ Gyr 
with the error bars derived from the variations in- 
troduced by the log{g) resolution. The radius of 



Eps Indi Bb, taken from Burrows et al. (1997), is 
~0.93 Rj. Generally, the theoretical radii do not 
vary by more than ~10%, with lower radii associ- 
ated with higher values of g. Figure[5] displays the 
measured spectrum together with the best fitting 
model. 

The fits to Eps Indi Ba are somewhat less good, 
but constraining. Fitting the Z and J bands was 
compromised by ambiguities in the cloud opacities 
and KI wings. To approximately achieve the ob- 
served flux level in the K band required clouds to 
suppress the fluxes in the Z and J bands. With- 
out cloud opacity, which naturally is most effica- 
cious in these bands, the Z and J fluxes would be 
too high, while the K band flux would simultane- 
ously be too low as shown in flgure[4j Glouds level 
out these bands and without clouds there are no 
good fits to Eps Indi Ba. An initial investigation 
showed that the best fit modal particle sizes were 
^^15 ± 5 /im. This is an interesting conclusion, 
albeit arrived at in lieu of a robust theory of con- 
densates in brown dwarf atmospheres. The three 
best fitting models are derived for (Tcs / log(g)) of 
(1250 / 5.3), (1250 / 5.2) and (1300 / 5.3) in order 
of preference corresponding to masses and ages of 
(55Mj / 1.7Gyr), (48Mj / 1.2Gyr) and (56Mj 
/ 1.5 Gyr). We finally estimate mass and age of 
Eps Indi Ba of 53+^^^ Mj and i.s+O-S -o.s 
The radius of Eps Indi Ba, taken from 'B urrows] 
et al.| (|1997|), is -0.83 Rj. Figure|4] displays the 



measured spectrum together with the best fitting 
model and a plot portraying a model, but with- 
out clouds, demonstrating the important role of 
clouds. 

There is a remarkable agreement in the best 
fitting ages of Eps Indi Ba and Bb of 1.5 Gyr, so 
the models are consistent with the assumption of 
coevality. However, the age range over which the 
models produce reasonable fits is of the order ±0.5 
Gyr, so in total the modelling results agree well 
with the age estimate of Eps Indi A of 



Gyr by [Lachaume et al., ( ,1999) suggesting overall 
spectral and evolutionary model consistency. Our 
mass estimates are somewhat higher than previous 
estimates (47 and 28 Mj for Eps Indi Ba and Bb, 
respectively, McCaughrean et al. (2004l), but still 
consistent with those considering the error bars. 

There is enough leeway in these fits that does 
not allow us to further constrain the ages and 
masses beyond this range. Ambiguities in con- 
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Fig. 4. — NIR spectrum of Eps Indi Ba. The un- 
derlying solid gray curve shows the best model fit. 
The dotted line shows a similar model not consid- 
ering clouds. 
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Fig. 5.— NIR spectrum of Eps Indi Bb. The 
underlying gray curve shows the best model fit. 



densate modeling and the age of Eps Indi A and 
degeneracies along the Tcs/g line in the spectral 
models are the culprits. Nevertheless, Eps Indi B 
is the most promising binary brown dwarf system 
with which to make absolute comparisons and for 
which astronomers can test the general theory of 
brown dwarfs. Hence, these data should be a goad 
to further measurements and a motivation for im- 
proved theory. 

4. Conclusions 

We have presented the 1st J-band as well as 
new H- and K-band spectra of Eps Indi Ba and 
Bb. By comparison with theoretical models we 
were able to derive strong constraints on effec- 
tive temperature, surface gravity and age of both 
brown dwarfs. The results confirm the approxi- 
mate coevality of the stellar and the two brown 
dwarf components of the Eps Indi system. The 
suppressed J-band and enhanced K-band flux of 
Eps Ind Ba indicates that a noticeable cloud layer 
is still present in a T1.5 dwarf like Eps Indi Ba 
at an age of ~1.5Gyr, whereas for an T6 dwarf 
like Eps Ind Bb clouds become much less impor- 
tant at the same age as clouds either sink below 
the photosphere (larger scale heights) or thin out 
(see, e.g., iCushing et al. (2008 1). Once dynam- 
ical mass estimates derived from monitoring of 
the orbital motion of the system become available, 
our slightly increased mass predictions (compared 
McCaughrean et al.| ( |2004) ) of SS+^-'^Mj and 



to 
34 



l^^Mj for Eps Indi Ba and Bb, respectively, will 
provide a crucial test for the validity of our mod- 
eling. Recent dynamical mass estimates point in- 
deed towards a significantly higher system mass 



for the Eps Ind Ba/Bb system (Cardoso et al. 
2008). Hence the Eps Indi system continues to 



be a benchmark object for testing and calibrating 
brown dwarf evolutionary and atmospheric mod- 
els. 
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